A new view of the spin echo diffusive diffraction on porous 

structures 

*Janez Stepisnik 
March 2002 



PACS: 76.60.Lz, 61.43.Gt, 66.30.Hs 



Abstract 



Analysis with the characteristic functional of stochastic motion is used for the gradient 
spin echo measurement of restricted motion to clarify details of the diffraction-like effect in a 
porous structure. It gives the diffusive diffraction as an interference of spin phase shifts due 
to the back-flow of spins bouncing at the boundaries, when mean displacement of scattered 
spins is equal to the spin phase grating prepared by applied magnetic field gradients. The 
diffraction patterns convey information about morphology of the surrounding media at times 
long enough that opposite boundaries are restricting displacements. The method explains 
the dependence of diffraction on the time and width of gradient pulses, as observed at the 
experiments and the simulations. It also enlightens the analysis of transport properties 
by the spin echo, particularly in systems, where the motion is restricted by structure or 
configuration. 
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1 Introduction 



Measurements of the molecular displacements through precession of their atomic nuclear spins 
in non- uniform magnetic field by the magnetic resonance spin echo[jl|] has gained a most decisive 
role at the studies of molecular transport within porous structures that may comprise system 
as diverse as sandstone rocks, catalysts, colloids, or biological tissue. The field has an extensive 
literature and a wide range of experiments has been performed. The methodology of the Pulsed 
Gradient Spin Echo (PGSE) [Q| has been successfully implemented to measure diffusion in sys- 
tems for which a constrained molecular motion caused a deviation from Fickian behavior. Based 
on the diffusion propagator formalism ||] , which gives the spin echo as a Fourier transform of the 
probability distribution , Callaghan has introduced the concept of diffusive diffraction of spin 
echo in a porous media pj. In this way, one can extract information not only about the motion 
but also about the morphology of the surrounding medium ^, pi 0]. Unlike NMR Imaging 
where molecular position are recorded to a resolution on the order of lOfim, the method is able 
to achieve the resolution of displacement measurement some two to three orders of magnitude 
better, and pushes the lower limit of NMR resolution into a nanometer range. 

However, the use of probability distribution function is just one of ways to consider the 
transport properties of spins in terms of the averaged characteristics of the elementary events. 
The method of characteristic functional, i.e. the Fourier transform of the probability distri- 
bution [10], is an alternative way to average the spin phase fluctuation. With the cumulant 
expansion in the Gaussian approximation, it gives the relation between the spin echo and the 
spectrum of the single particle velocity correlation function (VCF) JTl| , |l^ , 13 1. It provides a 
simple derivation of spin echo attenuation for any gradient pulse sequence. It also shows that 
a sufficiently fast and a properly shaped gradient sequence, like the Modulated Gradient Spin 



Echo (MGSE) |Tj, 15|, conveys information not only about macroscopic flow and diffusion, but 



also about the details of motion on the molecular level in the frequency range of gradient mod- 
ulation. The method has been successfully implemented to measure the diffusion spectrum and 
the flow dispersion in porous media [jl6L |l7| ]. 

This study is an attempt to apply this method to the analysis of the spin echo diffusive 
diffraction. It aims to clear up details like the diffraction dependences on time, the dependence 
on the width of gradient pulse and on the type of applied gradient sequences as obtained by 
experiments and computer simulations [18, |l9|, 20 1, but remains unanswered in the frame of 
propagator theory. The method might elucidate the dependence of diffraction on the mechanism 
of scattering at interfaces as well. 



2 Spin echo and average of spin phase fluctuation 

Whenever in NMR a non-uniform magnetic field is used to encode the spin magnetization for 
motion rather than position, the spin echo is used to refocus any spin phase shift, due to 
absolute spin position. Thus, the perturbations of spin phase, due to displacements r{t) of spins 
in the non-uniform magnetic field gradient G(t), can be written as 6(t) = 7 J Q T G(t) ■ r(t) dt = 
— Jq F(t) ■ v(t)dt, where F(t) = 7 / * G(t')dt' is a factor of spin dephasing, r is the time of 
refocusing, and v(i) is spin velocity. Regarding their location in a non-uniform magnetic field 
and rf exitation, one can distinguish subgroups of spins according to their precession frequency 



to obtain the spin echo as |21] 
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where Ej is the normalized amplitude and the the motion average of particles in the j-th sub- 
ensemble is (...) t. In the case of a sequence with two sharp gradient pulses of widths 5 and 
interspaced for A, i.e. a sharp PGSE sequence where F = q = 7#G, the mean of the spin 
phase fluctuation can be worked out with the conditional probability distribution (the diffusion 
propagator) P(r', t'\r, t), || |||, ||. It gives the spin echo as 

E(r,q) = ^^ (e^-( r i( T )- r i(°)) ^| P(r')dr , y'drP(r,r|r / ,A) e if l-( r - r/ ). (2) 
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Clearly, the spin echo appears as a Fourier transform of the probability distribution with respect 
to q Q. Thus, the measurement of E(r, q) gives an average propagator of restricted diffusion. 
This approach also explains the diffusive diffractions, which appear at long diffusion times and 
strong gradients ||. 

According to statistical physics, the stochastic process is determined knowing either the 



probability distribution or the characteristic function, which is the Fourier transform [22] of 
probability distribution, 

(e i f r ) = Je i f r P(r)dr, (3) 

The characteristic function exists even when the probability does not. By passing the stochastic 

properties from r to v, where r = / v(t')dt', the characteristic function is transformed into the 

Jo 

characteristic functional that can be worked out with the cumulant expansion method as 

i f f(t')v(t>) dt' i f f(t')( V (t')) df -\f ati f dt 2 m).(v(tiHPi))cf(t2) + ... 

(e Jo ), = e Jo 2 Jo Jo (4) 

where f(t) is an arbitrary function. With the cumulant expansion of characteristic functional, the 
stochastic process is determined by the expectation (v(t)), and by the set of correlation functions, 
of which the second order is {v(ti)v(t2)) c = { v (t) v (fy) ~ i v (t)) 2 - Gaussian approximation of the 
stochastic process is the expansion to the second term, which is often used to describe physical 
processes and can be justified in many processes in the magnetic resonance as well. 

These basic facts from the theory of stochastic processes leads to understanding that the 
average of the spin phase fluctuation can be resolved in different manners depending on the 
applied gradient sequence. With a sharp PGSE sequence, the position of spin-bearing particle 
r(t) is considered as a stochastic variable permitting to treat the spin phase fluctuation by using 
a probability distribution as shown in Eqj2|. In the cases of finite gradient pulses, multi-pulse 
gradient sequences or gradients of general waveform, the short gradient pulse approximation 
fails, and the stochastic properties can be passed to the velocity of spin, v. As the spin phase 
average has a form of characteristic functional, it can be treated by the well-developed methods 
of statistical physics. With the cumulant expansion in the Gaussian approximation, the spin 
echo appears as 

E( T ) = Y J E jo e i( t ) ^-f 3 ^ T \ (5) 
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where the phase shift depends on the mean local spin velocity 

^•(T)=-r pw^^)*, (6) 
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and the spin echo attenuation is related to the velocity correlation function 

&( r ) = \ f T r F ^* 1 ) • < v i(*!) v i&))c ■ Fj(t 2 ) dhdh. (7) 

^ JO JO 

For a free motion in a simple fluid, the mean velocity is assumed to be zero (v g j(t)) = 0, while 
VCF ca be approximated with (v g j(t) v g j(0)) = 2DS(t), due to a short memory of molecular 
collisions. Substitutions into Eq.|| gives a spin echo with exponential decay, which is in proportion 
to the mean squared displacement. According to Eq.||, its Fourier transform, with respect to 
parameter F(i) 2 dt, gives the probability distribution function for a free Fick's diffusion as 
expected. 

However, for motion in complex systems, there are a number of characteristic time-scales, 
which may correspond to frequencies in the time regime of spin echo. These include a long 
tail decay of VCF in liquids [p^, p4| , tube disengagement times in entangled polymers [25], a 
characteristic negative decay of VCF in confined fluids^, |2?J. Such times are more closely 
related to the structural dynamics of a liquid than to local particle motion. It usually appears 
as an anomalous time dependent attenuation of the spin echo. 

In the previous implementations of cumulant expansion method in the analysis of the spin 
echo measurement of restricted motions [^, |l^] , the cumulants were averaged over the volume of 
confinement. Therefore, the phase shift of spin echo has been neglected because of (v) = 0,and 
all attention was focused to the second term of expansion. It provides a known relation between 



the spin echo attenuation and VCF 17 ] . Such approximation is permitted as long as the gradient 
sequences are short enough or the gradients are weak enough that the spin displacement is short 
compare to the spin phase grating, ((F.Ar) 2 ) <C 1. Whenever it is not fulfilled, local distributions 
of the phase shifts and of the spin echo attenuation must be carefully considered. 

3 Distribution of velocity correlation function and mean veloc- 
ity of confined motion 

For the diffusion in restricted geometries, the approximations with (v g (t)) = and with VCF 
as a delta function are reasonable as long as the number of molecular impacts at walls is small 
compare to the number of intermolecular collisions, r c <C t <C t w . At longer times, we need 
better approximations. The solution of the Langevin equation for diffusion between parallel 
planes provides such VCF [27|, but an identical results one can obtain with the use of the 



probability distribution from the Fick's diffusion equation to average velocity correlations as 



shown in reference [29]. Therefore, we use this method to obtain the distribution of the mean 
velocity and VCF within the space of pore as follows. 

In the general form, the conditional probability function for a restricted diffusion is 

P(r, t\r ,t ) = Mr)Mro)e- k2 D ^ 

k 

, where k and eigen functions tpk( r ) characterize the compartment geometry. The mean spin 
velocity follows from the first derivative of the mean spin displacement as 

(g.v(r,t)) = j t J v g.(r-r')P(r,t\r')dr', (8) 

where the unit vector g is aligned along the magnetic field gradient. The second derivative of 
the mean squared displacement gives VCF as 

(g.v(r,t) v(r,t).g> = ^ J ( g .(r-r')) 2 P(r,t\r')dr' (9) 
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Figure 1: Distribution of the diffusive phase shifts between plan-parallel planes for the sharp 
PGSE 

Provided distributions make possible to obtain the spin echo for any gradient sequence by 
using relations of Eq|] and Eq.^. Here, we consider the simplest case, i.e. the application of 
the sharp PGSE sequence in order to enlighten the basic fact of diffusive diffractions from the 
point of view of new approach. With the eigen functions for the diffusion between plan parallel 
planes, Eqj6|, Eq]7| with EqJ^ and Eq.|| give the distributions for the phase shifts and the spin 
echo attenuation at different times as shown in Fig.[l]. At early times after the first gradient 
pulse, only spins in the proximity of wall are involved in the scattering. The component of 
mean velocity outward of boundaries gives the resulting phase shift that is in proportion to 
spin displacement, 0(r, r) ~ q\/Dt, but with the opposite sign at the facing planes. At longer 
time, when a number of scattered spins increases, the phase shift develops into almost linear 
dependence on location, 0(r, r) « q(r — (r)). As displayed in Fig.0, the distribution of the spin 



P(x,t) o 




x/d 



Figure 2: Distribution of the diffusion attenuation between plan-parallel planes for the sharp 
PGSE 

echo attenuation at early times is almost uniform but with small dales in the proximities of the 
walls. With time, it develops into a well-defined distribution with a maximum in the center of 
compartment. 

With such strong dependences on the location, the averaging over the space of compartment 
for the mean velocity and VCF is not permitted. It requires the analysis that includes the local 
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details of motion. Thus, the spin echo in the continuum limit has to be written 



E(t) 



E(r) e 



i 4>{t, r) - /3(r, r) dr 3 



(10) 



where 0(r, r) and (3(t, r) describe the distribution of spin phase and of attenuation, respectively. 

Time (Dx/d ) 



Log(E/E ) 




Gradient (qd/27i) 



Figure 3: Diffraction-like features of sharp PGSE as a function of time and gradient for diffusion 
between plan-parallel planes. 

Integration of Eq.[To|, for the case of diffusion between parallel planes and with a sharp PGSE 
sequence, gives diffraction-like patterns as shown in Fig.||. Diffractions exhibit a dependence on 
gradient magnitude q as well as on time r. At short times, the diffraction minima are shifted 
toward larger q, and depend on the spin displacement as 1\J Drq ~ n2-K. At displacements 
large enough that a spin is starting to experience scattering on the opposite boundaries, the 
minima appear at value of q a ~ n2 n where a is about the diameter of pore. It occurs at the 
displacements above 2Dt > 0.3a 2 according to Fig.||. At these times, the sharp PGSE obtains 
the form 

Jv 

where Mi = j v rdr and M2 = f v r 2 dr, and the diffractions can pass on an information about 
morphology. 

For the diffusion between parallel planes that are interspaced for a, it gives 



1 



E(t) 



,2sin(ga/2) 



(12) 



which becomes is in the limit of of small q 



E(t) 



(13) 



Although the last result is the same as that obtained with the propagator approach^], the 
dependence of the spin echo diffraction patterns on q as well as on the interval of measurement 
are very different as shown in Fig.|3[ 
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4 Conclusion 



The analysis with the cumulant expansion of the characteristic functional of stochastic motion 
provides some new details of diffraction-like features that occur at the gradient spin echo mea- 
surement of diffusion and flow within porous media. As in optics, the diffraction is interference 
of waves with different wave vector, here it is the interference of phase shifts because of back- 
flow of spins scattered at boundaries, whenever the repulsive displacements are comparable to 
the spin phase grating created by applied magnetic field gradients. At long times, when a spin 
starts to experience scattering on the opposite boundaries, the diffraction conveys information 
about morphology of the surrounding media. This approach is able to explain the diffraction 
dependence on the time and the duration of gradient pulses, as observed at experiments and 
simulations and can be implemented for any gradient sequence. It also casts a new light to 
the spin echo measurements of transport properties, particularly in complex systems where the 
motion is constrained by structure or configuration. 
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